Vortex and Droplet Engineering 

in 

§ Holographic Superconductors 

o 

(N 
< 

Tameem Albash, Clifford V. Johnson 

<N 

Department of Physics and Astronomy 
^ University of Southern California 

^ Los Angeles, CA 90089-0484, U.S.A. 

talbash, johnsonl, [at] usc.edu 

> 

in 

0^ 

^ Abstract 

^ We give a detailed account of the construction of non-trivial localized solutions in a 

2+1 dimensional model of superconductors using a 3+1 dimensional gravitational dual 
O theory of a black hole coupled to a scalar field. The solutions are found in the presence 

>■ of a background magnetic field. We use numerical and analytic techniques to solve the 

full Maxwell-scalar equations of motion in the background geometry, finding conden- 
sate droplet solutions, and vortex solutions possessing a conserved winding number. 
These solutions and their properties, which we uncover, help shed light on key features 
of the {B, T) phase diagram. 
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1 Introduction 



An holographic model of some of the key phenomenological attributes of superconductivity in 
2+1 dimensions was proposed in ref. [T]. It works roughly as follows (a more detailed review 
will follow in the next section). The dual is a simple model of gravity in four dimensions 
(with negative cosmological constant) coupled to a U{1) gauge field and a minimally coupled 
charged complex scalar ^. Asymptotic values of the scalar on the boundary correspond to 
the vacuum expectation value (vev) of a charged operator in the 2+1 dimensional theory. 
For high temperatures (relative to a scale set by non-zero charge density in the model) the 
system is in a normal phase represented by a charged black hole solution in the gravitational 
dual with the scalar set to zero. At a critical temperature the system undergoes a phase 
transition, the U{1) getting spontaneously broken by a non-zero vev of the charged operatoij^ 
The gravitational description of this is a charged black hole with a non-trivial scalar profile 
that gives the vev on the boundary. This means that the black hole has "scalar hair" in 
this regime. (For a discussion of violations of no-hair theorems in this context, see ref. [3]. 
There, it was shown that it is possible in this context, for large enough charges, equivalent 
to the low temperature regime here.) The authors of ref. [Tj showed using linear response 
theory that the DC conductivity of this new phase diverges in a manner consistent with the 
expectation that the system is in a superconducting phas^ The authors carried out further 
study of the system in ref. [2]. 

Our focus in this paper is the system in an external magnetic field, continuing the 
work we began in ref. [6]. Generically, for non-zero magnetic field B filling the two spatial 
dimensions, it is inconsistent to have non-trivial spatially independent solutions on the 
boundary, and we present and study two classes of localized solutions in some detail. The 
first is a "droplet" solution, the prototype of which was found in our earlier work[B| as a 
strip in 2D (straightforwardly generalized to circular symmetry in ref. [2]), and the second 
is a vortex solution, with integer winding number ^ G Z, which is entirely new. We obtain 
these as full solutions of the Maxwell-scalar sector in a limit, and determine a number of 
their properties. 

Our analysis in various connected limits shows where these solutions can exist in the 
{B,T) plane. There is a critical line below which droplets are not found, while vortices can 

^Strictly speaking, the U{1) that is broken is global on the boundary, but it can be gauged in a number 
of ways without affecting the conclusions. See e.g. ref. [2j. Note that a global U{1) does not restrict us to a 
spatially independent magnetic field. 

^Other holographic superconductors are available. See e.g. refs. [HIS]. 
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be found there. Our interpretation is that this region is the superconducting phase, and that 
for non-zero B, the vortices develop, trapping the magnetic flux into filaments, as is familiar 
in type II superconductors. Above the critical line, the system leaves the superconducting 
phase, and either forms droplets of condensate or simply reverts to the normal phase (dual 
to a dyonic black hole with zero scalar everywhere, which may well yield lower action than 
the droplets if we had back-reacting solutions to work with). 

In section 2 we review the model, and discuss the two limits in which most of our 
studies will be carried out. Section 3 reviews the spatially independent solution correspond- 
ing to the prototype superconducting solution. We briefly discuss our numerical approach 
to finding the solution as a warmup for the more difficult problems in the sequel. Section 4 
presents our search for and construction of the non-trivial spatially dependent solutions 
corresponding to condensate droplets and to vortices. We discuss the numerical methods 
we used to find them, and then examine a number of their properties. Section 5 examines 
aspects of the solutions' stability. We conclude in section 6, and we also present two ap- 
pendices. One appendix establishes the normalisation of our gauge/gravity dual dictionary, 
while the other discusses the flux quantization in our vortices. 

2 The Model, and Two Limits 
2.1 The Model 

The holographic model of superconductivity in 2+1 dimensions proposed in ref.[T] is a model 
of gravity in four dimensions coupled to a U{1) gauge field and a minimally coupled charged 
complex scalar \E' with potential V^(|\E'|) = — 2|\E'p/L^. There is a negative cosmological 
constant that defines a scale L via A = —3/L^. The action is: 

^buik = ^ / d'xV^S^R +^,+L' (^-\f' - \d^ - zgA^f _ y | , (i) 

where ^1 = SvtGn is the gravitational coupling and our signature is ( — h ++). 

We will use coordinates {t, z, r, (p) for much of our discussion, with t time, (r, 0) 
forming a plane, and z a "radial" coordinate for our asymptotically AdS4 spacetimes such 
that z = is the boundary at infinity. The AdS4 metric is: 

ds^ = ^ (dz^ - df + dr^ + r^d<f) . (2) 
z^ 

Note that the mass of the scalar m\ = —2/LF' is above the Breitenlohner-Freedman stability 
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bound m| 



BF ~ 



9/4L- 



^ for scalars in AdS4. We will write the scalar as: 



P 



(3) 



V2L 



Near the boundary z = we have: 



P 



Piz + P2Z 



,2 



(4) 



where pi [i = 1, 2) sets the vacuum expectation value (vev) of an operator Oi with dimension 
A = i [8J. Only one of these vevs can be non-zero at a time, and we will choose to study the 
case of i = 1, for much of the paper. Our charged operator will be the order parameter for 
the spontaneous breaking of the U{1) symmetry. A gauge field of the form A = Atdt does 
not give an electric field in the dual theory on (r, 0), but defines instead [9J a f/(l) charge 
density, p and its conjugate chemical potential p, as we will recall below. 

Black holes in this study will be planar, i.e., their horizons are an (r, 0) plane at 
some finite z = Zh- In the familiar manner, their Hawking temperature T and mass per unit 
horizon area e = M/V, corresponds to the dual 2+1 dimensional system at temperature T 
and with energy density e. Generically, the black hole will couple to the gauge sector, having 
some profile for the field At. The temperature T will have dependence on the charge density 
parameter p. This is quite natural since without the charge density there is no other scale in 
the theory, and there would be no meaning to a high or low temperature phase, and hence 
no possibility of a phase transition. 

The high T phase of the theory is simply the charged black hole (Reissner-Nordstrom 
(AdS-RN)) with the scalar \1/ vanishing. This corresponds to the non-superconducting or 
"normal" phase of the theory, where the order parameter vanishes. The mass of the scalar \1/ 
is set not just by but by the density p through the coupling to the gauge field. In fact 

m|, decreases with T until at it goes below m'^p, becoming tachyonic. The theory seeks a 
new solution, in which the black hole is no longer AdS-RN, but one that has a non-trivial 
profile for 

In studying the system in a magnetic field background, there are some generic expec- 
tations to consider. The magnetic field B (which fills the two dimensions of the supercon- 
ducting theory), also contributes to m|,, via its square, but contributes with opposite sign 
to the electric contribution of the background. It therefore lowers the temperature Tc at 
which m|, falls below m|p, triggering the phase transition. On these grounds alone one then 
(naively) expects a critical line in the {B,T) plane connecting (0,Tc) to some {Bc,0), but it 
is important to determine exactly what physics lies on either side of the line. Our solutions 
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and our study of their properties help in estabhshing some of this. The solutions can be 
found by solving equations of motion in certain limits and we remind the reader of them in 
the next two subsections. 



2.2 The Decoupling Limit 

After a field redefinition: 



^ , ^ ^ , (5) 



our action Q becomes: 



(6) 



If we consider the limit g ^ oo, then the Maxwell-scalar sector decouples from gravity. 
This allows us to work with a fixed uncharged background, which we take to be the AdS4- 
Schwarzschild (AdS-Sch) black hole, given by: 

ds' = ^ (-/ (z) dt' + dr' + rW) + KvT-^d^' ' (7) 
^ Z'^ J [z) 

where f {z) = 1 — z^. The coordinate z is a dimensionless parameter scaled such that the 
event horizon is at Zh = 1- The Hawking temperature is given by the usual Gibbons-Hawking 
calculus [To] : 



Note that a is related to the mass of the black hole: 



(9) 



where V is the volume of the (r, 0) plane. In terms of the two real fields (p, 6) into which we 
decomposed \1/ into in equation (|3| we have: 



-L' \d^ - + = -^G"^'^ [d,pd,p + p2 {d,9d,9 - 2A^d,9 + A,A,)] + j^p' . (10) 

From the action, we derive the equation of motion for the fields p, 6, and A^: 

[V^G'^'-d^p) - G^'^p {A^ - 0^6) {A, - d,e) + ^p = 0, 

-^d, (v^G^-p2 _ gj^^ ^ , (11) 

-^d, (v^G'^^G^-Fa.) - f^p2 (A, - d,9) = . 
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While the background itself will have no charge in this limit, there will of course still be a 
non-trivial gauge field A, and for an electric background A = Atdt, we will have, as z —>■ 0: 

— = At^fi-pz, (12) 
a 

defining a chemical potential and a charge density p. 
2.3 The Probe Limit 

Sometimes we will also work in a probe limit, where we take the scalar in the Maxwell-scalar 
sector to be small, and hence not back-reacting on either the geometry. In general, we can 
do this at arbitrary g. (We will combine this with the decoupling limit {g — > oo) for one 
case, as we shall see later.) For finite g we will consider our small non-backreacting scalar 
to be moving in a dyonic Reissner-Nordstrom background, given by [Hj : 

ds' = '^{-f{z)de + dr' + rW) + ^^,. (13) 

^ ' Z'^ J [z] 

F = 2ha'^rdr A dcp + 2qadz A dt , 

fiz) = 1 + {h' + g2) z^-(l + h^ + q^) z^ = {i-z){z^ + z + l- (/i^ + q^) z^) . 
The temperature and charge density are given by: 

^ = ^ = 4;^('-' ' p = wJW^) = -^'''' ■ ^'^^ 

We choose a gauge such that the gauge field is written as: 

A = ha^r'^dct) + 2qa {z - I) dt . (15) 

3 Spatially Independent Solution 

We begin by considering a spatially independent solution, reviewing the original presentation 



of ref. [1], working in the decoupling limit of section 2.2 We take an ansatz for the fields 
given by: 

9 = const , p = p{z) , At = aAt (z) , A^ = . (16) 
where p and At are dimensionless fields. The equations of motion are given by: 



d!p+[^-l)d,P + j,pA't + ^p = o, 



dlAt - ^p'At = . (17) 
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We can study the equations' behaviour near the event horizon {i.e. as z ^ 1): 



z=l 



A, 







2 = 1 



d^A, + ip^d^A 



. 



2 = 1 



We note that the only free variables (to be chosen) at the event horizon are p(l) (or dzp{V)) 
and dzAt{l). The other limit to study is to consider the behavior near the AdS boundary 
{z -> 0): 

1 



OzP - -Ozp + -^P 







J 2 = 



d'iAt - ^-p'At 



0, 



2=0 



which has as solutions 

p{z ^ Qi) ^ Piz + p2Z^ , At{z^Qi) 



p~ p z 



;i8) 



where pi, p2, p, and p are constants related to the vev of a A = 1 operator, the vev of 
a A = 2 operator, the chemical potential, and the charge density of the dual field theory 
respectively. The solution for p at the AdS boundary admits two normalisable modes, and 
therefore the constants are associated with vevs of two separate operators. Only one of these 
vevs is to be non-zero at a time, and the two different gauge theories are related to each 
other via a Legendre transformation [8]. 

3.1 Numerical Analysis 

To simplify the numerical analysis, it is convenient to define a new field R{z) such that: 

R{z) = zp{z) . (19) 



With this redefinition, the boundary condition of having either pi or p2 in equation (18) 



to be zero becomes the requirement of having either a Dirichlet or a Neumann boundary 
condition on R at the AdS boundary. The equations in the bulk of AdS are given by: 



dlR+i^dzR+j,RAl+(^f^ 



+ 



22/ 



/? = , 



dlAt - jR^At = 



(20) 



and we solve them using a shooting method (discretizing using finite differences) with shoot- 
ing conditions: 



R{1) = const , dzR{l) = -|i?(l) , d^Ml) = -^i?(l) - 3^ (9,A(1) ) i?(l) . 
Ml) = , 9,A(1) = const , d%il) = -\R{iydMi) ■ 



(21) 
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The solution at z = goes as: 

R{z ^0) = R1 + R2 z , At{z ^0) = fi- p z . (22) 

We fix -R(l) and then tune dzAt{l) until the solution satisfies the necessary Dirichlet or 
Neumann boundary condition at z = 0. We then read off the scalar and also the value 
of p for that solution, which defines the temperature. We can determine Tc since there is a 
minimum charge density (over temperature squared) needed for the scalar field to condense. 
Note that there are multiple choices for dzAt{l) that give the necessary boundary condition 
at the AdS boundary, sample solutions of which we present in figure [T] Solutions with a 
greater number of nodes are associated with higher chemical potential/charge density. These 
solutions are of a higher energy and so are thermodynamically unfavorable, therefore we only 
present results of the zero-node solutions in what follows. In figure |2] we show the solutions 

R(z) 

1.0 ^ 



0.5 



-0.5 



\ 0.2 9A _,0.6 0.8 1.0 



Figure 1: Three solutions with the same i?(l) but different dzAt{l) that satisfy Dirichlet boundary 
conditions at 2; = 0. The sokitions are distinguished by the number of nodes (times they cross the 
z-axis) they have. 

for the scalar values pi,2 at the boundary which give the vevs of the operators (9i,2- As 
anticipated, in each case, the vev of the operator is zero above T/Tc = 1. Below T/Tc = 1, 
it is not zero, showing the spontaneous breaking of the U{1) symmetry. 

4 Spatially Dependent Solutions 

A non-zero magnetic field B in the (r, 0) plane will correspond to some non-zero A^{r). In 
such a case, consistency of the solution requires the fields to have some spatial dependence 
in the plane. This situation was studied in the linear case in ref. [6], (see also ref. |2]) but 



here we consider the full non-linear problem of equation (11). First, notice that the U{1) 
gauge transformation acts as: 

p ^ p , e^d + X, A^^A,, + d^K. (23) 
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In the previous section we chose 9 to have no non-trivial dependence. Naively, it would seem 
that we can freely shift 6 by gauge transformations. However, this freedom is only available 
if the gauge symmetry is not broken. We will return to this once we have constructed the 
solutions. This motivates us to consider the following ansatz: 



= C + ^<P, P = pir,z) , At = aAt{r,z) , A^ = (r, z) 



(24) 



where we have defined a dimensionless radial coordinate f = ar, dimensionless fields p. At, 
and A^, and (CO ^'^^ constants where ^ is an integer. Under this ansatz, the equations of 
motion reduce to: 



a^i^ + i^d^A^ + ) {d',A^ - \d,A^ - (^^^ _ ^ , (25) 



dlAt^)[dlAt^\d,At 



^fP^At = , 



where the equation of motion for the field 9 is trivially satisfied by our ansatz. Near the 
event horizon, these equations reduce to the following conditions that must be satisfied: 



djA^ - \drA^ - [a^ - = 3a,i, 
it(f,2= 1) = , 



2 = 1 



dip = -p-lp{d,At 



2 = 1 



2 = 1 



d'A, 



A^-i)p^-2d,A 



2 = 1 



?,dlAt = ( djAt + \drAt - p'At ) + 2Atp^ 
9 



2 = 1 



(26) 



where in the first equation we have used that At{f, z — 1) — 0. We now have three free 
functions to fix in these equations, d^pir, z — 1), dzA^lf, z — 1), and dzAf{f, z — 1), which 
determine the spatial profile of the solutions at the event horizon. Note that in order to 
avoid a divergence in the equation for p at f = 0, we must have that near f — 0, the field p 
must go as r^. This motivates the following field redefinitions: 



p = zf^R{f, z) , = r^A , 



(27) 



where R near f = is a non-zero value. The particular redefinition of simplifies the 
numerical analysis. Our new equations of motion for these fields are: 

-2" ~2-^ ^""^^ "■^^^ 



d^^R + ^-d,R +\( dlR + '^^^dfR +^R-l;R (r^A - ' 



/ 



/ 



diA 



^-d,A+-^{dlA+^^dfA 



zf z^ ' z^f 



- ^ + 



R' A- 



R 



1 
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(28) 



dlR + '^^^dfR + - 4^ 



(f'A-i)- 



R = 3ai? 



2 = 1 



dlR ^-\r- 2dM - Ir (d.Aty 
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dU + ^drA - A 



z=l 



2=1 



2=1 



djf + -dff- r^^R^f = 6d,f 



2=1 



d^T = -2d,T 



2=1 



where we have used that: 



(29) 



lim At{r, z) = lim (1 — z) T{r, z) . 

2-+1 2— »1 



(30) 



In particular, at the event horizon, we can expand the fields near f = as: 

1 



lim i?(f, 2;) = Rq{z)[1 + -a2f^ + 0{f 



1. 



lim^(f,z) = Ao{z)\l + -h2r^ + 0{r)\ , 



limT (f, z) = To{z) ( 1 + -caf^ + ^(f^) 



(31) 
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Substitution into the equation of motion at the horizon yields: 

i?o(l) + 3a,i?o(l) - 2eAo(l)/2o(l) 



02 



2(e + i)i?o(i) 



C2 



Ao(l)flo(l)2+3a,Ao(l) 

4Ao(l) ' 
-gRo(l)2+3e,Ao(l) 

4Ao(l) ' 

4Ao(l) ' 
To(l)Ro(l)^+6a,To(l) 
2To(l) 
6azTo(l) 
2To(l) ' 



(32) 



The solutions we study are characterized by the value of ^ and the f asymptotic behavior 
of the field R. For any allowed value of ^, the solution for R can asymptote to zero or to a 
constant non-zero value. 

We will be extracting non-trivial profiles for the fields at the boundary at z = as 
follows: 



r A{r, z) 



fi{f) - p{f)z , 
acl,{r) + J^{r)z , 

p{f,z) = zR{f-,z) = R{f-,0)z + d^R{f,0)z^ , (33) 

where /z is related to the chemical potential, p is related to the charge density, J^j, is related 
to the azimuthal current, and is related to the magnetic field via = {df.a^)/f. For the 
exact relationships, please consult Appendix \K\ 

4.1 Droplet Solutions 

We consider the case of the case of = 0, and use Oi as our order parameter. For this choice 
of ^ the solution that asymptotes to a constant value is simply the spatially-independent 
solution described earlier in section [3j In this section we consider solutions that asymptote 
to zero. To that end, we fix the following functions to: 



--(1 + 7)%,!) 



9,i(f , 1) = 



(34) 



where 7 is a positive number. With these choices, the equations of motion at the event 
horizon reduce to: 



d^R + \dfR - r^A^R + 7^ 







z=l 



J 2=1 



d?A + ^dfA - R^A 



z=l 



d}.f + \dff - R^f = 6d,f 



z=l 



dlA 



z=l 



(35) 



2=1 
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The coefficients in equation (31) are given by: 



7 , Rojlf To(l)i?o(l)^ + 6g,To(l) 



4.1.1 Numerical Procedures 



We begin by solving equations (29), which are at the event horizon. For a given -Ro(l); we 
find that there is a specific value for Ao{l) and To(l) that gives regular solutions for the 
three functions R{r,l), A{f,l), T(f, 1). The coupled ordinary differential equations are 
discretized using an explicit finite difference method, and we determine the values of Ao(l) 
and To(l) using a shooting method. By this we mean that we pick values for Ao{l) and 
To(l) at the origin and "shoot" towards f — > fmax, where fmax is the largest radius out to 
which we will construct our solutions. Typically, this leads to a divergence in the functions, 
and therefore we iterate the procedure, fine-tuning our initial conditions such that a regular 
solution is found. This has now determined our initial conditions for the bulk problem. 

The initial conditions at the event horizon having been determined, we solve the bulk 



equations of motion (28) and shoot towards the boundary at z = 0. The coupled partial 
differential equations are discretized using a finite difference method, and we adjust the mesh 
spacings Az and Ar until we achieve stability for our code. 

In order to satisfy the necessary boundary conditions at the AdS boundary, we try 
to minimize the positive area under the curve of dzR{r, 0). We accomplish this minimization 
by fine tuning our choice of dzT{r, 1) at the event horizon. This is acheived by expanding it 
in an appropriate basis of functions in f and using a Monte-Carlo method to determine the 
coefficients (with the area playing the role of energy) . 

4.1.2 Sample Solutions 

To give a sense of how the solutions behave, we present multiple solutions for multiple values 
of 7 and -Ro(l)- The solutions are presented in figure [sj From figures 3(g) and 3(h) 



we see 



that for a given value of 7, for various initial condition values of the scalar field, we get the 
same asymptotic charge density. We learned from the spatially independent solution that 
the ratio of T/Tc is determined by the charge density, and therefore we learn that 7 fixes the 
value of T/Tc- In fact, as 7 — 0, we have T/T^ 1 and as 7 — > 00, we have T/T^ — > 0. From 
figures 3(c)| and |3(d) , we see that the magnetic field asymptotes to a constant value, which 



indicates that the solutions "live" in a background magnetic field. As the magnitude of the 
scalar field increases, the value of this background magnetic field rises. Note also that how 
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(e) Current density for 7 = 0.5 
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(g) Charge density for 7 = 0.5 
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Figure 3: Droplet solutions for 7 = 0.5 on the left and 7 
T/Tc ^ 0.84 and T/T^ ^ 0.67 respectively. 



2 on the right. They correspond to 



the magnetic field behaves in the core of the droplet varies considerably between low and 



high temperatures. At low temperatures (figure 3(d)), the magnetic field is enhanced by the 
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droplet, whereas as high temperatures (figure 3(c) ), the droplet weakens the magnetic field in 
the "core." This variation in behavior suggests that these are perhaps not superconducting 
droplets as was thought (Gj [2]. 

We can try to study the minimum magnetic field needed to first form these droplet 
solutions. This would correspond to studying the problem in the limit where the magnitude 
of the scalar field is approaching zero, i.e. the perturbative or probe limit. In this limit, it 
is consistent to take A to be a constant in both f and z and to take At to only depend on z. 
The equation of motion for R with A = 7/2 reduces to: 



dm + -dfR - -f'YR + iR 

r 4 

This equation has (lowest energy) solution given by: 

R{f,z = 1) = i?o(l)exp 

Next, we can solve for the z dependence by solving: 

dlR + i^d^R - )^R + j,RAl + 

dlAt = , 



. 



(37) 



z=l 



-7- 



f2 



(38) 



i? = , 



(39) 



with the appropriate boundary conditions at the AdS boundary. This in turn fixes the value 
of the temperature of the solution. The value of the magnetic field found here corresponds 
to the critical magnetic field at which the droplet solutions first form. We draw the corre- 
sponding diagram in figure |4j The source of the divergence in the critical magnetic field as 




Figure 4: The Hmiting droplet line in the g ^ 00 limit. Below this, droplets disappear. 

T/Tc ^ is clear. As the magnetic field and charge density grow in value, they begin to 
back-react on the geometry, and our decoupling breaks down. However, our results allow 
us to shed new light on the calculation done in ref . [6] . There, the prototype solution in the 
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droplet class was first uncovered in the probe limit. Now we see that the probe limit is the 
correct limit to study the onset of the droplets. So far in this section we've seen them in the 
decoupling g ^ oo limit, and further taking the probe limit gives a limiting line at which 
they drop to zero height, ceasing to exist for lower B. 

We can study this further (and extend to lower T /T^) by working again in the probe 



limit, but at arbitrary as outlined in section 2^ Here, the background is now a charged 
black hole solution, our method of taking into account some of the back-reaction of the 
gauge fields. In this limit the p equation of motion is given by: 

\2 o 



f 



1 



(i-zY 



^9'q'p + —~p = ^. (40) 



This equation has a separable solution that we write as: 



p = zZ{z)R{r) 



(41) 



and we write their respective equations of motion as: 



d^R + ^dfR - ^ (V/i^) R - 2ghR = 



diZ + '—d,Z - -2ghZ + 



f f 
The solution for R{r) is given by: 



:i-^) 



Ag\^Z + 



2 f 2 

z^f fz 



i?(f) = exp {-ghf^/2) 



Z = 



(42) 
(43) 

(44) 



We can solve the equation for Z{z) using the same arguments as before, and we get the 
solutions shown in figure |5} Our claim is that the perturbative scalar field on the dyonic 
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Figure 5: Limiting droplet line for three values of g with T = (3 — q"^) /Air. 
black hole describes the entire phase transition line in the phase diagram, and the problem 
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studied earlier is simply the g ^ oo limit near T/Tc — 1. To prove this, we do several non- 
trivial checks. First, we check that the dyonic theory can predict the critical temperature of 
the g ^ oo theory. In order to do this, we make the identification: 



= Pg- 



(45) 



by comparing the charge density of both theories. Next, we know that as the magnetic field 
approaches zero, the droplets appear at the critical temperature (or T/Tc = 1) for both 
theories. Given that we defined Tc in different ways for both theories, by setting them equal, 
we should be able to calculate the relationship between and ^/p that we saw in the g ^ oo 
theory. In particular: 

(46) 



1 



Tc/a = — {3-ql) = o^^gq^ . 

Therefore, we can solve for a g ^ oo. We present the results in figure |6j As one can 
see, in the limit of (7 — > oo we indeed recover the values 0.226 and 0.118 respectively, which 
were obtained earlier in the g ^ 00 probe case (see caption of figure pi). We can also check 
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Figure 6: Convergence to the g ^ 00 critical temperature. The curves asymptote to 0.225492 and 
0.118412 respectively. See text for discussion. 

whether the phase diagrams coincide. This is presented in figure [7j As we see, the dyonic 
black hole results very quickly approach our results for the g ^ 00 case for a range of T/T^ 
near one. We are now in a position to answer what happens when T/Tc — >■ for the g ^ 00 
limit. The zero temperature limit requires us to take: 



q^ + h^ = 3. 



(47) 



As 5^ — > 00, we find that regularity of any solution at zero temperature requires us to take: 



3 1 
2g 



33/4 



1 



(48) 



Therefore, the Gaussian profile in equation (44) vanishes in the limit of zero temperature 
and (7 — > 00 (note that this does not happen at finite g). Therefore, the droplet no longer 
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Figure 7: The droplet limiting curves for a range of couplings, after rescaling to include the g ^ oo 
case. Tc = aa^2qg. 



exists in that Umit. Another way to see this is that in this Hmit, the dimensionless quantity 
gBzjTl as used in figure [t] diverges. 

4.2 Vortex Solutions 

We now consider solutions for the scalar that asymptote to a constant non-zero value. These 
are the vortex solutions. 



4.2.1 Numerical Procedures 

The numerical procedure for the vortex is almost identical to that of the droplet (see section 



4.1.1). We found it much more difficult to solve for the initial functions on the horizon using 
a shooting method, and so we inserted an initial guess function for an approximation to the 
scalar field -R(r, 1) at the event horizon, parameteriszed by two constants -Rq and -Ri. We 
then use that function to solve for the field A(f^\\ however it turns that for a given v4o(l) 
and -Ro) there is a specific -Ri that leads to a regular solution; we again use a shooting method 
to determine this constant. The constant To(l) is determined in the same way. With both 
i?(f, 1) and y4(f, 1) determined, we have fully fixed 9^-R(f, 1). The bulk shooting problem is 
tackled as before. Note that the leading term in the expansion of d^T is the constant already 



determined in the spatially independent problem using equations (21), since our vortices 
asymptote to that case. 

4.2.2 Sample Solutions 

Here, we illustrate the case of ^ = 1 and ^ = 2 and again use 0\ as our order parameter. We 



again consider the equations of motion given in equations (28) and (29). For simplicity, we 
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focus on the case of d^A = 0. We note that as the solutions approach a constant value, they 
should asymptote to the spatially-independent solutions we have presented earlier. This 
in turn allows us to define the temperature at which a given solution exists. We present 
examples of such solutions in figures [8] and [9j In figure |8| we see that very far away from 
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Figure 8: Vortex Solutions for ^ = 1 (LHS) and ^ = 2 (RHS). 



the origin, the scalar field has a constant vev, but as it approaches the core of the vortex, it 
decreases to zero value. The behavior at the origin is of course determined by the choice of ^. 

and 



The choice of ^ also influences how the gauge field behaves. In figures 9(g' 



9(h) 



we 



see that the value of A^ asymptotes to the value of ^. This is exactly the behavior required 
for the magnetic flux penetrating the vortex to be quantized with value 27r^. We review 
this briefly in Appendix [B} Indeed, ^ defines a non-trivial topological winding number: The 
scalar p becomes constant at infinity, breaking the U{1). Therefore the gauge symmetry of 
equation (23) cannot be used to unwind 6. Gauge symmetry is unbroken at infinity for the 



droplets, so ^ is not a winding number for them. The current density J(j){r) (figures 9(e) and 



9(f) ) is zero asymptotically and peaks in a ring around the core, supporting the magnetic 



field, as expected for a vortex. 



In figure |9(b) , we find that the charge density near the origin begins to oscillate as 
the density drifts downwards. It is difficult to say whether or not this is a physical attribute 
of the solutions or whether it is an artifact of our scheme to find the appropriate shooting 
functions that satisfy the z = boundary condition. If they are physical, they may be caused 
by screening effects being strong in the core of the vortex. It is interesting to note that this 
behavior appears to be absent for the ^ = 1 vortex (see figure 9(a) ), although we do see that 
there is a transition from the charge density increasing in the core to decreasing in the core 
as the temperature is lowered. Another curiosity for the ^ = 1 vortex is the behavior of the 



magnetic field near the origin (see figure 9(c) Instead of fiattening out as is the case for the 
^ = 2 solutions, it dips slightly downwards. We expect this also to be a numerical artifact. 
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Figure 9: Vortex Solutions for ^ = 1 (LHS) and ^ = 2 (RHS). 
since the ^ = 1 case is more numerically sensitive because of its sharper profile near f = 0. 
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5 Stability and Quasinormal Modes 

In this section we consider the stabihty of our solutions to fluctuations in the fields. We work 
in the decoupling limit and proceed by considering fluctuations about our classical solutions: 

^^ = 4°)+A«, p = pW + p«, ^ = + (49) 

where the flelds with superscripts (0) are the background flelds and the fields with super- 
scipts (1) are the fiuctuation fields. Note in particular that because we are working in the 
g —>■ oo limit, we do not consider fiuctuations of the metric since those come at 0{l/g^). 
The action to quadratic order in the fluctuations is given by: 

s - d'xV^ {-Ig^' + (aw - - d^e^^y) 

+4p(°)p(i) (4°) - (4) - dj^'^) + (pW)' (4) - d,9^'^) (4) - dj^'^) 

+TA~p'''T-TKl^n^) ■ (50) 

The resulting equations of motion are given by: 

(y^ (2p(°)p(i)G'''' (4°^ - + (pW)'G'''^ (4^ - a.^(i)))) = o , 

5.1 Spatially Independent Solution 

We can consider the ansatz: 

AW^ae-'-'A^^\oj,z). (51) 
The equation of motion for the fluctuation fleld is given by: 

dlA^^' + jd.A'^' + 1^4^ - y {r''') ' = . (52) 
Near the event horizon, the fleld satisfles an equation of the form: 

S'A" - Y^/A'' + = • (53) 
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This has solutions given by ingoing (negative) and outgoing (positive) waves: 

i« oc (1 - zf''^ . (54) 

The appropriate condition to have at the event horizon is of ingoing waves, and therefore, 
we make the following field redefinition: 



iW = (l-z 
The equation of motion is now given by: 

dlx 



(55) 



1 



+- 



lUJ 



-dzX 



1 1 

72^+^--, ZX 



. (56) 



Expanding near the event horizon, we get the following restrictions on the initial conditions: 



dzX 



2 "1 



9 — Qiu) 



-X 



(57) 



J 2 = 1 



By requiring that x{z = 0) = 0, we find that this condition is only satisfied for discrete 
values of uj. We present some of the values in table [l} In particular, since the imaginary 
part of UJ is negative, this corresponds to having fluctuations that decay away. Therefore, 
the constant solutions are stable under fluctuations. 



n 


uo 





1.42804 - 1.87689 i 


1 


3.07885 - 2.79650 i 



Table 1: Values of uj that give regularizable solutions at T/Tc = 0.772. 



5.1.1 Hydrodynamic Limit and Conductivity 

As an aside, we can push our analysis a little more to compute the DC conductivity (already 
done in ref. [1]). We proceed by studying the problem in the limit where uj approaches zero. 
We therefore consider solving for the field xi^^ z) as an expansion in uj: 

x{.^,z) = xo{.z)+uJXi{.z) + ... (58) 
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As explained in ref. [12], in the hydrodynamic limit and where the field is solved as an 
expansion in cD, the normalizibility condition cannot be satisfied by the terms Xi- We proceed 
by focusing on the xo term. The equation of motion is given by: 

d!xo + jd,Xo-j{R^'^) Xo = 0. (59) 
Note that at the AdS {z = 0) boundary the solution for xo is given by: 

limxo=( ^^ 'tJt' (60) 

where and are constants and is proportional to the current. In particular, if we 
define the conductivity using: 

^ Jx Jx (61) 

Ex iuiAx 

Therefore for T < Tc, Im(cr) oc co~^, and therefore by the Kramers-Kronig relations, we have 
that Re(cr) oc 6{uj) as shown in ref. [T]. 

5.2 Spatially Dependent Solutions 

We begin by slightly simplifying by considering the following field behavior: 

p(i) ^ e-'^^p^'^\u, r, z) , aJ^^ = ae-''^^A[^\u, r, z) , A^^^ = e-''^^A^^\u, r, z) , (62) 
A^^^ = ae-'''^Ai^\u,f,z) , Ai^^ = ae-'^*A^}\Co,f,z) , 9^^^ = . (63) 

where we have defined dimensionless fields and variables uj = au. This reduces the equations 
to: 



-l(p(°))^i(^) = o 



-^(p(°))^i« = 0, (64) 
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and 

-;^(P'"')''41" = 0. 



+ (if)' + -^ifil^'p'"' + ^p'" = . (65) 

Consider two possible ansatze. First, let us consider the radial gauge choice A^P = 0. In 
this gauge, the equation for gives the following restriction: 

(a^i^) + ^i(^)) = -JdJ^^ . (66) 

This result is interesting because it suggests that in the limit oi u ^ 0, there is a consistent 
solution with Ar^"* = 0. This suggests that for this ansatz, Ai'^^ oc uj. However, implementing 



the restriction in equation (66) is not trivial, and therefore, we consider a different ansatz. 
For lack of a better name, we call the temporal gauge, A[^'' = 0. In this gauge, the equation 
for ij^^ gives the following restriction: 



+ J (d.A^^ + U^') ) = JyP«P^°)iS°^ . (67) 



In particular, for this ansatz, in the limit of — > 0, we see that we must have p^^-* oc u which 
in turn means that we must have A^^^ oc uj for consistency. This particular restriction is 
more straightforward to implement since we can directly insert it into the equations for A^z'^ 
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and Ai^\ The resulting equations of motion are given by: 

-P7(^^°^)'4'^=0' (68) 

and 

What is of particular interest is that we find that the equations for p^^^ and A^^'^ only depend 
on each other, and so we can solve for these two fields independently of A^r^ and A^P . For 
simplicity, we therefore set these fields to zero. We now note that under this assumption, the 
equation of motion for A^P appears to be completely independent of A\^\ In particular, one 
could imagine first solving the equation for A^ and then substituting the solution into the 
equation for A^ and solving for A^ under the constraint described above. Unfortunately, we 
do not know how to solve the partial differential equation with the constraint, and therefore 
we choose to solve only for the field A^ . To proceed, we would like to note that near the 
event horizon, all the fields satisfy the same equation: 

dlx + + = . (70) 

This equation has solutions given by in-going (negative) and out-going (positive) waves: 

X oc (1 - ;2)^*"/' . (71) 

The correct boundary condition to have at the event horizon is of in-going waves, so we 
redefine the field to reflect this: 

iW = (l-^)-^/^x(r,^) . (72) 
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The equation of motion for Z is now given by: 
1 /' 



Six 



9 fl 



z) J 



dzX + -pX + 



lUJ 

y 



l-z 



){dh + \d,x-{R^'^)'^ =0 



At the event horizon, this gives: 



(73) 



r 



X + iuJX + ^uj x + (2iu; - 3) d^x 
1 



^^X + ( 2 - -i^ ) d,x + (a; - 6z) X = 



(74) 



The first point we need to address when solving these equations is what the appropriate 
choice for oj is. Since for the droplet solutions the scalar field asymptotes to zero, we expect 
that in order to satisfy the normalizibility conditions for large f, the correct choice for u are 
the quasi-normal solutions when there is no scalar field present. We present some of the 
values in table [2j Using these values, we try to find normalizable solutions for the droplet 



n 


UJ 





0.66252 - 2.60058 i 


1 


2.20850 - 2.51647 i 


2 


4.03511 - 2.45694 i 


3 


6.06406 - 2.41322 i 


4 


8.26149 - 2.37694 i 



Table 2: Values of u that give regularizable solutions for zero scalar field. 



scalar field solutions. What is interesting is that we can only find solutions for ^{uj) > 4 
(note that this does not mean that solutions for smaller ^{uj) do not exist, just that our 
extensive numerical search could not find them). We present an example of this solution 
in figure 10 Therefore, these results suggest (note that this result is not conclusive since 
we have not solved for the field Ai^'') that the droplet solutions are stable under quadratic 
fluctuations. The situation for the vortex solutions is very similar, except that now, the u) 
values we use are those of the spatially independent solutions described earlier. Once again, 
we find that solutions can be easily found above the lower lying modes, although we find 



them to be much higher above as illustrated by the example shown in figure 10(b) 
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(b) Linear Vortex Fluctuation with 
w = 13.12 - 2.32i and T jT^ = 0.77. 



Figure 10: Solutions for the fluctuation fleld x- 



6 Conclusion 

We have constructed two broad famihes of locahzed solution to the equations of motion for 
ref. [l]'s holographic model of a superconductor, and considered several of their key properties. 

The vortices, with winding number ^, contain 27r^ units of magnetic flux, and are 
candidates to fill out the superconducting part of the phase diagram in the presence of an 
external magnetic field. This is because a lattice of them in the (r, 0) plane can trap flux lines 
of an applied external S-field into filaments as it passes through the two dimensional sample. 
Vortices presumably repel each other, and so such a lattice will cost energy. Therefore at 
some critical Bj(r) the system will seek a lower energy phase, possibly returning to the 
normal phase. We have not constructed such a lattice, and further study to understand such 
a configuration is a very interesting avenue of research to pursue. Forming such a lattice is 
a method by which, at a given T/Tc < 1, the superconducting phase can be made to persist 
in some constant background 5, even though the system cannot eject the magnetic field 
entirely a la the Meissner effect (there is a nice energetic argument in ref. [2] as to why the 
Meissner effect is not possible in this two dimensional case). This vortex phase is of course 
the same method by which a standard type II superconductor can persist beyond the (lower) 
critical line at which the Meissner effect disappears, and we expect that it applies here. The 
study of multi-vortex solutions needed to establish this is left for further study. 

Crucially, we've established that the droplets do not exist below a certain critical value 
of B, dropping to zero height on a family of lines that we were able to compute explicitly. For 
this and a variety of other stated reasons, and also considering the fact that they are of finite 
size and hence a lattice or gas of them would not give a connected superconducting path for 
charge transport, we believe that they do not represent a superconducting phase. (Hence, we 
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disagree with the statements made about the phase diagram in ref. [2]. The authors find the 
critical fine, but state (similarly to our ref. 0) that the droplets exist below the line, and are 
superconducting. As they did not have the full droplet solutions, nor the vortex solutions, 
their analyses are not sufficient to make these determinations.). They seem to represent a 
non-superconducting phase that is inhomogeneous. Whether or not the droplets are the 
favoured solution for arbitrarily large B is an interesting question. There is the possibility 
that the system may prefer to return to the normal phase represented by a dyonic black 
hole with zero scalar. Our partial stability analysis showed that the droplets we studied are 
stable against ffuctuations of the fields, but there is the possibility that other fiuctuation 
modes may be undamped. We mention here that we also noticed the curious fact that the 
droplet solution at higher magnetic field (that we presented earlier) contains regions where 
the local value of the squared scalar mass is below the Breitenlohner-Freedman bound. 

It would also be of interest to establish whether the critical line where the vortex phase 
would disappear coincides with the hmiting line where the droplets' existence begins. We 
conjecture this to be likely on the grounds that we have found no other candidate solutions 
to fill an intermediate region. While this is the simplest possibility, further study is needed 
to establish it firmly. 
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A Normalizations in the Holographic Dictionary 

We recall the AdS dictionary (working in Euclidean metric): 

(^exp j (j)oO^ = exp (-S'on-sheii[0o]) = Z . (75) 

Taking derivatives on both sides with respect to the boundary source 0o gives us our AdS 
dictionary: 

(0(xi) . . . 0(a:„)) = iPV)-- }^^,^''jJr-, ■ ■ ■ jJr-^^ ■ (76) 

00^0 d(po{Xi) d(f)o{Xn) 

We define the the free energy density of the dual theory to be given: 

^ = ^yS'on-shell ■ (77) 



27 



where (3 is the inverse temperature and V is the "spatial volume" of the dual theory. If we 
use the notation that we are working in AdS^+i, then V has mass dimension d — 1. For us, 
d = ?>. In the Euclidean language, our action is given by: 



(7J 



where we emphasize that although the metric is now purely positive, we are still using At 
and now a Wick rotated version of it. In the dual theory, the charge density is given by: 

where /i is the chemical potential and x represents the space-time coordinates in the dual 
field theory. Using the AdS/CFT dictionary, we can write: 

5J-' 1 (J^on-shell 



6fi{x) (3V6At{x,0) 



(80) 



Using the action given in equation (78), we find: 



g (^Son-shell "^^ ^ Q /i / n\ ^'^ ^ 2 a A ^ n\ 



pVSAt{x,0) 24 g ^ ' ^ 24 g 
where we have used that: 

|M£;i^,w<«. ,82) 

and we have dropped the contribution coming from the event horizon. Therefore, the end 
result is given by: 

1 

p{x) = --^-a^dMx,0) . (83) 
Zk^ g 

A similar procedure allows us to calculate the vev of the azimuthal current as well: 

Note that here we are using the vector field that has been rescaled by g, hence why the factor 
of g appears at the beginning of our definition. Next, we can calculate the form of the vevs 
of the A = 1 and A = 2 operators. To proceed, we take the A = 1 operator to be the source 
of the A = 2 operator [8]. Therefore, we write: 
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where we are using the notation that: 

z ^ 0) = zpi{x) + z^p2{x) . (86) 

To proceed, we calculate the variation of the bulk action (keeping only divergent and finite 
terms) : 

55'buik = -lim [ d^x L^a^ \-pi{x)5pi{x) + 2p2{x)5pi{x) + pi{x)5p2{x) \ . (87) 



The first term in parentheses is divergent, but it is removed by an appropriate counterterm: 



ScT = -:z-^^^iTO[i—— dx-p{x,z) , (88) 



which leaves us with: 

-2 



^^buik + SScT = -TT^^ / d^x L^a^p2{x)6pi{x) . (89) 
Zk,^ g J 



Therefore, we find: 



- ^J^) = 2,47 ''^"^ = ''^"^ ■ 

To calculate the vev of the A = 1 operator, we need to perform a Legendre transform on JF 

g = -^-^^ff- 1 Pi(^)P2(x) , (91) 

and we now have: 

Therefore, in order to satisfy the conditions: 

= -jr^frT\ ^ ((^i(^)) = -7z|^, (93) 

o{Oi[x)) b{02[x)) 



we choose: 



(Oi(x)) = --i—Lapi{x) = —i—api{x) , (94) 
y/ 2gK4 \/ 2gK4 

{02{x)) = -J^La'p2{x) = -^a'p2{x) . (95) 
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B Flux Quantization 



Let us review the flux quantization condition. Deep in the superconductor, away from the 
vortex, our resuhs indicate that we have = ^. This can be written as: 

^A^ = ^d^9 . (96) 

It is useful at this point to reahze that A^p/r transforms exactly as the vector A. Therefore, 



it is convenient to write equation (96) as: 

Ve = A. (97) 
Consider drawing a circle around the vortex, deep in the superconductor. We can choose to 



integrate equation (96) along the circle: 



V9-de = j A-M . (98) 

where A is the corresponding 1-form. The left hand side gives 27r^, and using Stokes' theorem 
on the RHS, we have: 

2ni = j {v X A^ da = j da {^drA^ . (99) 

Note that the quantity in brackets on the right is exactly the magnetic fleld, so we find 
that the flux of a single vortex is quantized to 27r^. 
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